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1.0  Summary 

The  development  of  a  room  temperature  (<50C)  fuel  cell  that  could  use  a  generally  available 
fuel  such  as  JP8  would  be  most  valuable.  However  there  are  no  known  catalysts  that  can 
selectively  activate  the  Carbon-carbon  (CC)  bonds  of  such  fuels  at  such  temperature.  The  goal  of 
this  project  was  to  investigate  whether  it  is  plausible  to  develop  such  a  catalysts.  To  do  this  we 
used  validated  first-principles  (quantum  mechanics)  based  simulations  to  investigate  the  barriers 
for  alkyl  carbon-carbon  cleavage  for  two  prototypical  systems: 

•  A  metal  alloy  catalyst 

•  An  organometallic  cluster  catalyst; 

that  might  serve  as  an  anode  for  electrochemical  power  generation. 

For  the  metal  alloy  catalyst,  we  used  a  Top-down  approach  where  we  determined  the  bond 
energies  of  the  alloy  catalyst  to  various  fuel  fragments  necessary  to  achieve  room  temperature 
decomposition  of  butane  and  oxidation  to  C02. 

For  organometallic  cluster  catalyst  we  used  a  bottom-up  strategy,  detennining  specific  metals, 
ligand  compositions,  solvents  and  operating  conditions  to  evaluate  which  conditions  could  lead 
to  reasonable  rates  at  room  temperature. 

For  the  bottom  up  approach,  we  used  the  density  functional  theory  (DFT)  form  of  quantum 
mechanics  to  investigate  homogeneous  catalysts.  We  decided  that  Rhodium  (Rh)  is  the  best 
metal  for  a  low  temperature  oxidation  of  Carbon-hydrogen  (C-H)  and  C-C  bonds.  We  then 
considered  a  number  of  possible  mechanisms  and  reaction  pathways  for  butane  oxidation.  The 
lowest  barriers  we  determined  for  these  mechanisms  are  30  kcal/mol,  which  is  probably  too  high 
for  reasonable  reaction  rates  at  room  temperature  (probably  requiring  -200C).  We  then  began  a 
study  of  Rh/CF  complexes,  such  as  spin  singlet  RhC^Cfi)2'  and  the  spin  doublet  Rh111  superoxide 
complex  RhCl3(H20)(02)1", which  we  consider  would  be  accessible  under  reaction  conditions. 

We  conclude  that  these  might  serve  as  3-  or  4-electron  oxidants.  Further  research  is  needed  to 
draw  a  firm  conclusion  about  the  prospect  of  these  and  other  possible  metal  complexes  to 
provide  suitable  catalysts. 

For  the  Top  Down  studies  we  used  ReaxFF,  the  bond  order  based  reactive  force  field 
developed  by  our  group  at  Caltech  in  reactive  molecular  dynamics  simulations  of  hydrocarbon 
decomposition  on  metal  surfaces.  This  was  used  to  model  the  hydrocarbon  fuel  conversion  on 
metal  surfaces,  where  instead  of  modeling  a  specific  system,  say  FeNi,  we  detennined  the 
surface-adsorb  ate  bond  energies  required  to  operate  at  room  temperature.  These  surface- 
absorbate  bond  energies  can  be  used  to  scan  the  periodic  table  to  find  which  combinations  of 
metals  is  most  suitable  for  successful  alkane  oxidation  at  low  temperature.  These  results  show 
that  no  pure  metals  are  suitable  as  catalysts  for  the  C-C  bond  cleavage  at  room  temperature.  This 
analysis  suggests  metal  alloy  systems  that  might  satisfy  the  desired  parameters  that  could  be 
investigated  further. 

Our  conclusion  is  that  there  is  some  possibility  of  developing  a  low  temperature  fuel  cell 
capable  of  using  butane  as  a  fuel.  However,  some  compoents  of  JP8  would  not  be  reactive  under 
these  conditions.  We  have  organized  a  conference  at  Caltech  to  criticize  this  analysis  and  to 
suggest  alternatives.  This  involves  six  speakers  on  the  morning  of  August  9  and  13  all  day  on 
August  12.  In  addition  there  will  be  3  additional  experts  at  both  sessions  to  provide  discussion. 
The  program  is  attached  in  section  4. 
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2.0  Introduction 

For  complete  electrochemical  oxidation  of  alkanes,  the  prospective  catalytic  systems  should 
display  as  many  as  possible  of  the  following  features: 

•  A  low  barrier  (probably  <15  kcal/mol)  for  activating  C-H  bonds  of  saturated  alkenes  to 
form  new  M-CFI3  and  M’-H  bonds 

•  A  low  barrier  (probably  <15  kcal/mol)  for  activating  C-C’  bonds  of  saturated  alkenes  to 
form  new  M-C  and  M’-C’  bonds 

•  A  mechanism  and  a  low  barrier  (probably  <15  kcal/mol)  for  forming  C-0  bonds 

•  Stability  against  oxidation  and  protonation  at  a  pH  and  potential  corresponding  to  a  low 
overpotential  above  the  effective  fueWCCF  oxidation  potential 

•  Not  poisoned  by  CO 

In  addition  electro/thennochemical  relationships  should  include 

•  Alkane-activating  catalyst  (“M-X”)  not  oxidized  at  operating  potential 
but 

•  Catalytic  intermediates  (many  electrons  worth)  that  all  can  be  oxidized 
and 

•  Catalytic  intermediates  (many  protons  worth)  acidic  at  the  operating  pH. 

(see  Figure  1) 


Figure  1.  Electrochemistry  features  of  a  catalyst  for  alkane  oxidation. 
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We  applied  top-down  and  bottom-up  strategies  to  evaluate  compositions  and 
structure-performance  relationships  that  might  satisfy  these  constraints. 


3.0  QM  Studies  of  Homogeneous  Catalysts  for  Oxidation  of  C-H  and  C-C  Bonds 


3.1  Methods,  Assumptions  &  Procedures 


In  order  to  identify  catalysts  relevant  for  the  tasks  of  the  project,  we  investigated 
computationally  many  possible  candidates  over  a  wide  range  of  pH’s  (with  roughly  correlated 
oxidation  potentials).  Free  energy  surfaces  for  complete  or  partial  catalytic  cycles  featuring  CH 
or  CC  cleavage  were  computed  using  density  functional  theory  (DFT)  coupled  to  a  polarizable 
continuum  solvation  model. 

We  took  a  homogeneous  (bpy)2Run(OH)2  complex  as  a  potential  catalyst  that  features 

•  a  later  transition  metal  with  the  ability  to  fonn  strong  M-C  bonds  and 

•  oxidative  stability  in  basic  aqueous  solution  above  the  CH4/CO2  oxidation  potential. 

To  avoid  oxidation  of  (bpy)2Ruu(H20)(0H)  to  (bpy)2Rum(OH)2,  the  complex  must  operate 
below  -0.4V  vs.  Standard  hydrogen  electrode  (SHE)  at  pH  =  14  (Rum  is  incompetent  for  H3C-H 
activation).  We  find  that  AG  for  CH4  activation  is  32.4  kcal/mol,  and  13.7  kcal/mol  uphill 
(Figure  2).  This  would  provide  reasonable  rates  at  200C.  Unfortunately,  the  calculated  barrier  for 
CH4  activation  (32  kcal/mol)  was  too  high  to  permit  ambient-temperature  CH  activation. 
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Figure  2.  The  he  free  energy  surface  for  stable  intermediates  as  a  function  of  pH  and 

electrochemical  potential. 


A  promising  system,  tungsten  phosphine  complex,  was  reported  by  Parkin  et  al.  [1]  to  cleave 
an  aromatic  C-C  bond  in  quinoxaline  stoichiometrically. 
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M-CH,  oxidation 


Figure  3.  Low-valent  organometallic  tungsten. 

To  determine 

1 .  whether  this  reaction  would  apply  to  general  hydrocarbon  substrates  or  is  specific  to  just 
the  quinoxaline  substrate  and 

2.  if  general  to  hydrocarbons,  whether  the  complex  is  stable  with  respect  to  oxidation  and 
compatible  with  C-0  bond  formation, 

we  calculated  the  free  energy  surfaces  for  the  interaction  of  the  (PMeijsW0  complex  with  both 
quinoxaline  and  butane.  This  allows  us  to  explore  the  scope  of  this  reaction.  Indeed  we  found 
that  the  reaction  of  this  highly  reducing  tungsten  complex  with  quinoxaline  is  8  kcal/mol 
exergonic,  but  unfortunately  the  similar  insertion  of  the  tungsten  into  a  C-C  bond  of  butane  is  70 
kcal/mol  endergonic  (Figure  3). 

This  suggests  that  the  reaction  is  not  a  generally  viable  mechanism  for  CC  cleavage.  A 
second  problem  with  this  catalyst  is  that  regenerating  a  metal  center  as  reducing  as  W°  is  likely 
quite  unfavorable  under  oxidizing  conditions. 
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As  a  third  class  of  systems  we  considered  the  catalytic  mixture  of  RI1CI3,  Cf,  I',  CO  and  O2 
in  a  water/fluorobutyric  acid  solvent.  This  system  was  reported  by  Sen  [2]  who  observed  that  it 
oxidizes  light  hydrocarbons  at  80-90°C.  This  system  displays  several  attractive  characteristics: 

•  low  temperature  oxidation  of  hydrocarbons, 

•  tolerance  of  CO,  and 

•  a  preference  for  CC  cleavage  over  CH  cleavage  that  is  most  unusual  in  higher  hydrocarbons. 
Drawbacks  in  the  Sen  system  include 

•  dependence  on  CO, 

•  the  price  of  rhodium,  and 

•  the  dependence  on  O2,  which  cannot  be  used  as  a  stoichiometric  oxidant  in  a  fuel  cell  setting. 
We  first  examined  the  RI1CI3  system.  Our  initial  computations  regarding  this  system  showed  that 

•  Rh11  chloride  salts  are  stable  in  a  mixture  of  Rh1  and  Rhm,  so  that  metallo-radicals  may  be 
present  and  have  a  role  the  catalytic  cycle, 

•  Oxidation  of  the  Rh1  chloride  to  Rh11  is  exothermic  even  at  SHE,  and 

•  Moreover,  oxidation  of  Rh11  to  Rh111  is  also  exothermic  (reflecting  the  driving  force  to  fonn 
Rh-Cl  bonds). 
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Figure  4.  AG  for  some  possible  reactions 
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These  free  energies  suggest  Rh1  and  Rhm  will  comproportionate  to  two  Rhn’s.  This  would  be  the 
resting  state  for  metallo-radicals  at  some  low  potentials,  suggesting  that  this  system  could  be 
useful  for  attacking  strong  CH  and  CC  bonds.  However  in  the  presence  of  O2,  Rh111  is  the 
thennodynamic  resting  state. 

We  calculate  that: 

•  C-H  cleavage  in  butane  by  insertion  of  Rh1  chloride  into  a  CH  bond  is  1 1.0  kcal/mol 
endergonic  (uphill),  with  a  26.3  kcal/mol  barrier. 

•  C-C  cleavage  by  insertion  of  Rh1  is  9.3  kcal/mol  endergonic  (uphill),  but  with  a  52. 1 
kcal/mol  barrier. 

•  C-H  cleavage  by  Rh111  is  16.5  kcal/mol  endergonic  (uphill)  (Figure  4). 

•  Nonradical  insertion  of  Rh1  into  the  C2-C3  bond  of  butane  yielded  a  calculated  barrier  of 
52  kcal/mol  showing  that  this  is  not  a  viable  pathway. 

However  this  very  complex  system  has  a  large  number  of  potential  reactions,  many  of  which  we 
have  not  yet  studied.  Thus  the  real  barriers  might  be  lower, 

Since  the  RhCl3/02/C0/I-/H20/F7C400H  catalyst  system  described  by  Sen  [2]  displays 
selective  methane  oxidation  and  butane  oxidation  to  C2  products  at  moderate  (90°C) 
temperatures,  we  selected  this  system  for  our  further  investigations.  Our  goal  was  to  detennine 
the  mechanism  by  which  C-C  bonds  are  broken  and  which  of  these  reaction  steps  directly 
involve  O2  and  whether  this  might  be  reproduced  in  an  anodic  setting.  Non-radical  mechanisms 
explored  earlier  had  barriers  too  high  to  be  responsible  for  reaction.  Thus  we  conclude  that  the 
critical  reactions  are  concerted.  There  are  many  combinations  of  Rh(III)/02/iodine  that  could 
serve  as  a  strong  oxidant  and  several  permutations  of  electron  transfer/hydrogen  abstraction/C-C 
cleavage  reactions.  Since  there  are  a  large  number  of  combinations,  investigating  them 
exhaustively  would  be  a  major  undertaking  too  extensive  for  this  project.  Indeed  despite  the  very 
promising  results  there  has  not  been  a  successful  experiment  reported  that  gives  any  suggestion 
of  the  critical  catalytic  steps. 

Therefore  we  first,  determined  the  energies  of  RI1/O2/I  complexes  and  their  1-,  2-,  3-electron 
reduction  potentials.  Then,  using  the  complexes  that  are  the  strongest  oxidants,  we  determined 
transition  states  for  C-C  cleavage.  Finally,  we  determined  whether  the  other  required  steps 
(electron  transfer,  hydrogen  atom  transfer)  are  energetically  favorable. 
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One  possibility  is  that  C-C  cleavage  in  butane  observed  by  Sen  [2]  involves  electron  transfer 
from  alkane  to  an  activated  oxidizing  complex.  To  detennine  what  magnitude  of  oxidation 
potentials  could  be  created  by  combining  rhodium,  dioxygen  and  iodide,  we  modeled  RI1/O2/I- 
complexes  with  different  bonding  arrangements  (Figure  5).  We  then  added  electrons  and  protons 
to  these  to  predict  what  magnitude  of  oxidation  potentials  could  be  generated. 

We  took  RhIIICl3l(C0)(H20)1'  as  a  representative  starting  point.  The  initial  coordination  of 
dioxygen  to  Rhm  in  this  case  (and  other  cases  we  have  studied)  is  unfavorable  in  free  energy  but 
required  for  activation.  O  bonds  are  present  in  both  the  dioxygen  complex  and  products, 
suggesting  a  role  for  iodine  in  the  O2  activation  process.  The  most  stable  of  the  three  isomers 
involves  stabilization  of  singlet-02  in  a  Rh-O-O-I  peroxide  fragment  (“G  =  25.2  kcal/mol).  After 
reduction  by  one  electron,  the  lowest-energy  state  we  were  able  to  find  that  keeps  the  0-0  bond 
intact  is  G  =  6.7  kcal/mol  endergonic  from  O2.  This  is  consistent  with  an  oxidation  potential  of 
0.8  V  vs.  SITE  for  the  reduction  of  the  Rh-O-  0-1  fragment  at  pH  =  0.  If  some  bond 
rearrangement  is  allowed,  a  more  stable  hypoiodite  complex  can  be  accessed  after  1-  electron 
reduction,  which  is  G  =  1 1.9  kcal/mol  downhill  from  O2.  While  it  is  feasible  that  the  superoxide 
complex  (G  =  6.7  kcal/mol)  can  rearrange  to  this  hypoiodite  complex  (G  =  -1 1.9  kcal/mol),  it  is 
not  obvious  that  there  is  a  favorable  pathway.  It  could  be  that  the  system  gets  “stuck”  at  the 
superoxide  state  prior  to  further  reduction.  After  the  first  reduction  has  generated  the  hydrogen 
superoxide  complex  (“G  =  6.7  kcal/mol),  the  next  reduction  is  much  more  favorable.  Breaking 
the  0-0  bond  (and  trapping  a  hydroxyl  radical  with  iodide)  upon  reduction  generates  the 
hypoiodous  acid  complex  (“G  =  -46.7  kcal/mol).  The  reduction  from  G  =  6.7  kcal/mol  to  G  =  - 
46.7  kcal/mol  implies  an  oxidation  potential  of  2.3  V  vs.  SHE.  This  is  formidable,  and  well 
higher  than  the  overall  O2  +  2e"  +  2H+  <->  2H2O  potential  of  1 .23  V.  Since  1-0  bonds  are  present 
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in  both  the  dioxygen  complex  and  products,  this  suggests  a  role  for  iodine  in  the  O2  activation 
process. 

Next,  we  determined  which  C-C  cleavage  mechanisms  are  available  from  Rhin/C>2/I  mixture 
and  found  that  all  located  pathways  have  inaccessible  barriers  (Figure  6). 


cu 


co 

'N 

ci 


G  =  0.0  kcal/mol 


cr 


°c 


CO  OH  -1 
CI3Rh"  '\ 

No^° 

H 

G  =  34.0 


CO 

a3Rhl“'^l..-o 

G  =  28.4 


e*  , 


co 
1  ■ 

CbRlV1-  n — 1 
XT  147 
H 

G  =  22.6 


-2 


'OH 


CO  " 

I  ^-1 

CI3Rh-  O 

0  1.42 
G  =  16.0(18.7) 


CO 

I 


-1 


CI3Rhl"''1°y! 


G  =  25.2 


Direct  C- 


ch3 

U  2 

co  Hj 

I  m"OH2 
CI3Rh'"  ,.-0 

G  =  82.2 


-1  r 


co 

«io  -CC 
.  G  =  'm  X2H3 


age  TS’s 


"c- 

h2 


CO 

ciXi"'‘OH 

.h‘°" 


G  =  395.1 


CO 


-1 


CI3Rh 


h2 


-OH, 


\ 


1  ""S;  h2 


o 

G  ~  '50 


CO 

Cl3RhC^OH 

OT34 
G  =  6.7 


CO 

I  ,ro— I 

CI3Rh 

'"OH, 

G  =  -26.9 


?-C  cleava 
H-ato 


vage  ll  s  f 
om  abstrac 


following 
action 


H,C 


CO 

CI3RhIM' 

h2"< 
H2C— -cl 


-1 1 


-OH 


'CH2 
G  =  378.5 


CO 

I  _^OH* 

CkRh1" 

-'O 

1N\c^c-”0 
.  . h2  H2 

G  ~  256 
CO  2 

I  mi-oh2 
CkRh1" 

No-"0 — I 


CO 


I 


-OH2 


clhRN0--OH 

H2  G  =  x 


G  =  *15.2 


H2C-..., 


04 


CO 

H3c.  CkRh^"1 
H2C— -c^" 

CH, 


-OH 


G=253.0 


Figure  6.  C-C  cleavage  mechanisms  available  from  RhHI/02/I  mixture. 


For  this  reason,  we  studied  a  similar  process  of  O2  coordination  and  reduction  using  a  Rh1 
starting  point.  The  electron-withdrawing  CO  present  stabilizes  the  Rh(I)  oxidation  state.  The  2- 
electron  reduction  of  O2  is  more  than  twice  as  favorable  as  the  one-electron  reduction  of  O2 
(Figure  7).  Unlike  Rhm,  Rh1  itself  can  provide  electrons  to  the  O2  fragment  (in  addition  to 
electrons  supplied  by  the  hydrocarbon.).  Reducing  the  Rh'-CU  complexes  by  one  electron  is  more 
exergonic  than  reducing  the  Rhm-02  complexes  by  one  electron  because  the  Rh1  center 
simultaneously  donates  a  second  electron  to  O2  to  make  the  total  reaction  a  two-electron 
reduction  of  CF. 
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Figure  7.  Thermo/electrochemistry  of  the  oxidant  RhVcVI  mixture. 


3.2  Results  &  Discussion 

Overall,  we  are  looking  for  low-barrier  processes  for  forming  alkyl  peroxide  intermediates, 
following  an  initial  CH  cleavage.  We  consider  that  a  RhV02  complex  or  a  Rhn/02  complex  may 
be  responsible  for  the  CH  cleavage,  since  they  showed  the  largest  reduction  potentials  in  the 
previous  studies.  Since  none  of  the  CH  cleavage  transition  states  using  the  RhVCF  are  consistent 
with  such  a  low-temperature  reaction.  The  best  hypothesis  for  CH  cleavage  involves  a  radical 
generated  by  Rh  and  O2  (Figure  8a)  (since  free  radicals  are  not  consistent  with  experimental 
observations).  Since  C-C  cleavage  barriers  beginning  from  Rh111  are  high,  we  consider  that  an 
alkyl  peroxide  may  be  formed  instead.  The  lowest  C-C  cleavage  barrier  we  have  been  able  to 
find  is  the  Acid-catalyzed  Baeyer-Villiger  decomposition  from  the  alkyl  peroxide  (Figure  8b). 
However,  we  must  show  that  the  alkyl  peroxide  can  be  formed  from  O2  +  Rh.  Electrochemical 
C-C  cleavage  is  possible,  but  since  this  involves  O2  directly,  it  is  not  useful.  Schulz  found  Co111 
ions  in  acetic  acid  oxidized  n-butane,  suggesting  methyl  ethyl  ketone  as  an  intennediate,  which  is 
also  consistent  with  a  BY  decomposition  [4]. 
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We  found  that  the  Rh1  complex  (CO^Rh^TFAXH^O)  (TFA=trifluoroacetate)  can  react  with 
butane  via  activation  of  a  secondary  CH  bond  with  a  thermally  accessible  barrier  of  24  kcal/mol. 
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Figure  8.  Possible  processes  involving  the  RhVcVI  mixture. 

Oxidation  of  the  resulting  Rh1  alkyl  complex  using  O2  is  exergonic,  and  we  identified  a 
nonradical  mechanism  for  C-C  cleavage  on  the  resulting  Rh111  alkyl  complex.  However, 
depending  on  the  ligand  set,  we  found  that  the  C-C  cleavage  reaction  leads  to  barriers  of  30-40 
kcal/mol,  showing  that  this  is  not  the  mechanism  of  the  Sen  catalyst.  This  mechanism,  in  which 
an  ethyl  group  migrates  1 ,  2  across  the  C-Rh  bond,  would  generate  a  Rh  ethyl  ethylidene 
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fragment.  This  could  subsequently  undergo  thermally  accessible  nucleophilic  attack  by  oxygen 
nucleophiles,  resulting  in  the  desired  C-0  bonds. 

Comparing  the  RhI/m  redox  couple  using  different  ligand  sets,  we  found  a  role  for  carbon 
monoxide.  The  most  stable  Rh1  complexes  contain  two  CO  ligands  and  are  the  only  Rh1 
complexes  whose  oxidation  to  Rh111  by  Vi  O2  is  endergonic.  Since  CO  is  known  to  poison  many 
catalysts  by  binding  to  active  sites,  we  sought  other  electron-withdrawing  ligands  that  could  be 
bound  to  Rh  statically  to  stabilize  the  low  oxidation  state  without  inhibiting  the  reaction.  We 
found  that  Fluorinated  phosphines  are  sufficiently  Ti-accepting  to  satisfy  this  role. 

In  our  next  step,  we  wanted  to  determine  the  mechanism  by  which  C-C  bonds  are  broken  in 
the  Sen  system.  Non-radical  mechanisms  explored  earlier  had  barriers  too  high  to  be  responsible 
for  reaction. 

Since  many  combinations  of  Rhm/02/iodine  could  serve  as  a  strong  oxidant  and  several 
permutations  of  electron  transfer/  hydrogen  abstraction/C-C  cleavage  reactions,  we  first 
determined  the  energies  of  RI1/O2/I  complexes  and  their  1-,  2-,  3-electron  reduction  potentials. 
Second,  using  the  complexes  that  were  the  strongest  oxidants,  we  looked  for  transition  states  for 
C-C  cleavage.  To  date,  however,  the  search  for  reactions  which  couple  CH  or  C-C  cleavage  to 
reduction  of  the  Rh/Cb/I-  fragment  has  only  yielded  reactions  with  thermally  inaccessible 
barriers. 

Unlike  Rh111,  Rh1  itself  can  provide  electrons  to  the  O2  fragment  (in  addition  to  electrons 
supplied  by  the  hydrocarbon.)  Reducing  the  RhI-02  complexes  by  one  electron  is  more  exergonic 
than  reducing  the  Rhin-C>2  complexes  by  one  electron  because  the  RhI  center  simultaneously 
donates  a  second  electron  to  O2  to  make  the  total  reaction  a  two-electron  reduction  of  O2.  The 
initial  reaction  of  butane  with  the  reactive  state  of  Rh  could  result  in  direct  CC  cleavage  or  CH 
activation  followed  by  C-C  cleavage  (Figure  9).  Rebound  of  an  alkyl  radical  resulting  from  CH 
cleavage  (G  =  -6  kcal/mol)  will  likely  form  a  Rh111  alkyl  (G  =  -27  kcal/mol),  since  this  rebound 
will  have  a  much  lower  barrier  (if  any)  than  ethylene  liberation  (25  kcal/mol).  New  transition 
states  for  C-C  cleavage  via  rearrangement  of  the  Rhm-alkyl  were  located  (G  =  12.2  and  24.5 
kcal/mol),  but  these  pose  inaccessible  barriers  of  40  kcal/mol  or  higher. 
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Figure  9.  C-C  cleavage  barriers  involving  the  RhVc^/I  mixture. 
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Our  calculations  showed  that  if  a  secondary  alkyl  peroxide  can  be  fonned,  a  Bayer- Villiger 
migration  can  occur  with  a  26  kcal/mol  activation  barrier,  that  may  be  sufficiently  low.  In  this 
process  the  C2-C3  bond  is  broken  as  an  ethyl  group  migrates  1,2  from  carbon  to  the  peroxide 
oxygen.  However  a  mechanism  for  fonning  the  isobutyl  peroxide  must  be  identified, 
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Our  computations  explored  possible  reactions  between  RhVCk/I-complexes  and  butane. 
Hydrogen-atom  transfer  from  butane  to  an  as-yet  unidentified  reactive  complex  of  oxygen  could 
yield  a  butyl  radical  and  a  Rh11  complex  exergonically.  The  likely  fate  of  such  a  pair  is  simply 
recombination  to  fonn  a  Rh111  alkyl  complex,  without  further  bond  cleavage.  We  did  find  a 
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competing  transition  state  for  C-C  cleavage  in  which  a  primary  butyl  radical  is  attacked  at  the  C3 
position  by  the  Rh11  metallo-radical,  resulting  in  C2-C3  cleavage  and  release  of  ethylene  and  a 
Rh111  ethyl  complex  (Figure  9).  The  path  poses  an  accessible  barrier  of  ~25  kcal/mol. 

We  also  investigated  the  activation  of  a  CH  bond  in  methane  by  Rh11  complexes 
(CO)Rhn(TFA)2Cl-  and  found  that  Rhu(TFA)2Cl-  yielded  inaccessible  activation  barriers  of  37 
kcal/mol  which  is  too  high  to  be  observed  thermally  (Figure  10).  The  TOF  for  methane  oxidation 
suggests  a  barrier  ~  28  kcal/mol.  These  barriers  are  similar  to  barriers  predicted  for  Rh(III)-TFA 
complexes. 

To  rationalize  the  predominance  of  C-C  cleavage  products  over  CH  cleavage  products  in  the 
oxidation  of  butane  by  Sen’s  catalyst,  Chepaikin  et  al.  [5]  proposed  that  CH  cleavage  occurs 
first.  But  the  resulting  catalyst  fragment  “X”  is  so  electrophilic  that  it  immediately  abstracts  an 
electron  from  the  alcohol  highest  occupied  molecular  orbital  (HOMO),  a  C-C  a-bond.  The 
resulting  radical  cation  “C+-C”  is  susceptible  to  nucleophilic  or  oxidative  attack  resulting  in  C-C 
cleavage.  This  requires  a  3-  or  4-electron  oxidant  “X-0”(Figure  1 1).  Thus,  we  propose  that  the 
initial  reaction  of  a  RhXX  complex  with  butane  involves  a  one-  or  two-electron  reduction  to 
result  in  another  reactive  complex  with  sufficient  oxidation  potential  to  abstract  an  electron  from 
the  central  C-C  bond  of  the  butyl  fragment  present.  The  oxidation  leads  to  C-C  cleavage.  This  is 
consistent  with  the  observations  that  O2  (a  4-electron  oxidant)  is  required  for  catalysis,  while 
traditionally  “more  reactive”  oxidants  (H2O2,  mCPBA,  IO4',  C6H5IO)  are  incompetent.  The 
ionization  potential  of  butanol  (calculated  2.4  V  vs.  SHE)  is  lower  than  that  of  butane  (calculated 
3.5  V).  Oxidation  of  butane  and  2-butanol  leads  to  a  large  weakening  and  lengthening  of  the  C-C 
bond.  Oxidation  of  the  alkane  (3.5  V)  would  be  exergonic,  but  if  the  oxidant  and  alkane  are  in 
proximity,  electron  transfer  coupled  to  bond  rearrangement  may  be  possible  (Figure  12). 
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Proceeding  from  the  above  discussion,  we  began  to  seek  RI1/O2  complexes  that  would  be 
accessible  under  reaction  conditions  and  which  could  potentially  serve  as  3-  or  4-electron 
oxidants.  We  identified  One  candidate  as  the  singlet  rhodium  chlorides,  RhCl3(02)  ’,  with  a 
singlet  wavefunction  similar  to  that  in  ozone  (Figure  13).  Its  formation  from  O2  and  Rh1  is 
exergonic,  and  its  wavefunction  suggests  it  may  be  a  competent  oxygen-atom  transfer  reagent. 
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We  assume  here  that  the  Rh(I)  reactant  did  not  “cost”  any  energy  to  fonn.  Since  Rh(I)  is 
observed  in  solution  in  the  presence  of  the  electron-withdrawing  CO,  this  is  likely. 
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Figure  13.  New  oxidized  Rh  chemistry. 


One  of  the  4  available  reductions  of  O?  is  “used  up”  when  O2  bonds  to  Rh11,  but  3  units 
remain,  and  the  radical  character  on  the  superoxide  ligand  is  likely  kinetically  reactive  (Figure 
13b).  We  assume  that  the  Rh11  reactant  here  did  not  “cost”  any  energy  to  form. 

Neither  the  Rhv  nor  RhIV  is  formed  exergonically  from  O2  (Figure  13c).  These  could  be 
formed  using  the  stronger  oxidant  H2O2  generated  in  situ.  However,  neither  of  these  are  3-  or  4- 
electron  oxidants. 

Another  candidate  might  be  the  doublet  Rh111  superoxide  complex  RhClXfUOXCX)1".  The 
radical  character  of  the  superoxide  unit  may  impart  reactivity  toward  hydrogen-atom  transfer 
reactions.  The  higher  oxidation  state  complexes  Rf^ClXOX^O)1"  and  Rh^OX^OXOH)1'  are 
less  stable  than  Rh111  precursors. 

Ultimately,  Rh-CX  complexes  like  those  in  Figure  13a  and  b  are  the  most  likely  active 
species,  being  formed  exergonically  from  the  reagents  present  and  capable  of  multi-electron 
oxidations. 


3.3  Conclusions 

In  summary,  we  find  that  the  “Sen”  system  [RI1CI3/I-/O2/CO  in  F3CC(0)0H/H20]  is  unique 
in  oxidizing  light  alkanes  at  85°C.  C-C  and  C-H  bonds  are  broken  and  C-0  bonds  are  formed.  In 
the  reaction  of  butane,  C-C  cleavage  is  faster  than  C-H  cleavage. 

Our  interest  is  to  identify  the  role  of  O2  and  CO  in  this  mechanism.  Since  they  will  not  be  present 
at  an  anode  they  must  be  replaced  by  viable  surrogates. 

Phosphine  ligands  stabilize  Rh1  (relative  to  Rh111)  so  that  it  may  continue  to  activate  C-H 
bonds  under  oxidizing  conditions.  Thus  phosphines  can  be  considered  as  a  replacement  for 
carbon  monoxide.  There  is  no  evidence  that  O2  can  be  replaced  for  forming  H3C-O  bonds:  the 
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methanol  produced  contains  only  labeled  oxygen  if  labeled  O2  is  employed,  and  calculated 
barriers  are  high,  ~30  kcal/mol.  However,  it  can  be  replaced  for  forming  bonds  between  oxygen 
and  heavier  alkyl  groups:  barriers  for  nucleophilic  attack  of  anions  on  RhIU-C2H5  are  below  20 
kcal/mol. 

4.0  Top-Down  determination  of  Metal-adatoms  bond  parameters  suitable  for  low 
temperature  fuel  cells  using  butane  as  the  fuel. 


4.1  Methods,  Assumptions  &  Procedures 

The  ReaxFF  force  field  is  a  first  principles  based  force  field  trained  to  describe  energetics 
and  reaction  barriers  to  match  an  extensive  set  of  quantum  mechanics  reaction  pathways.  To 
consider  the  parameters  that  a  suitable  metallic  alloy  must  have  to  serve  as  the  anode  for  a  low 
temperature  fuel  cell,  we  used  ReaxFF  to  model  the  hydrocarbon  fuel  conversion  on  metal 
surfaces  in  a  top-down  approach.  We  started  with  the  ReaxFF  already  developed  to  describe  the 
anode  reactions  of  hydrocarbons  for  Ni/SOFC  catalyst.9'11,12'14  Then  we  tuned  the  M-H,  M-C, 
and  M-0  parameters  to  reduce  the  reaction  temperature  to  -100C.  Then  we  used  these 
parameters  to  suggest  which  elements  to  consider  in  the  alloy.  Note  that  the  high  temperatures 
are  used  to  achieve  reaction  rates  sufficiently  high  that  we  could  get  results  in  a  day  or  so.  More 
realistic  temperatures  would  take  longer  and  would  be  used  for  the  final  alloys  selected. 

First,  we  used  ReaxFF  to  describe  reactions  of  butane  on  a  Ni  (1 1 1)  single  crystal  surface. 
Here  20  butane  molecules  were  exposed  to  a  Ni  slab  containing  480  atoms  (Figure  14).  The 
temperature  was  increased  from  50K  to  1250K  over  25ps  and  held  at  1250K  (NVT)  for  lOOps. 
The  cleavage  of  C-H  bonds  (with  Ni-C  bond  fonnation)  was  observed  after  87  ps,  but  C-C  bond 
cleavage  was  not  observed. 


Snapshot  at  87.5  ps 


Ni:  black; 
C:  red: 

H:  green 


Figure  14.  Butane  on  the  Ni(lll)  surface. 


Then  we  continuously  heated  the  system  from  1250  K  to  1500  K  over  100  ps  and  found  two 
cases  of  C-C  bond  breaking:  the  first  one  was  observed  at  1415  K  and  the  other  at  1500  K 
(Figure  15). 
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Figure  15.  Reaction  pathway  for  C-C  bond  breaking  at  T  =  1415  K. 


To  determine  how  strong  a  Ni-C  bond  should  be  in  order  that  a  surface-bound  butyl  fragment 
would  undergo  CC  cleavage  spontaneously,  we  modified  the  parameters  of  the  ReaxFF  force 
field.  The  Ni-C  single  bond  energy  was  increased  until  a  butyl  fragment  bound  to  the  Ni  surface 
through  Cl  and  C2  decomposed  without  a  barrier  into  methyne  and  propyl  fragments.  A  60% 
increase  in  the  Ni-C  bond  energy  was  required.  This  adjustment  of  the  ReaxFF  resulted  in  the 
effective  increase  of  the  Ni-CFfi  single  bond  from  43  to  85  kcal/mol,  and  in  the  increase  of  the 
Ni-CH  bond  (3  Ni-C  single  bonds)  from  123  to  257  kcal/mol.  After  that,  the  Ni*  metal  with  a 
bond  energy  to  a  methyl  fragment  42  kcal/mol  stronger  than  that  of  real  Ni  could  break  a 
dehydrogenated  C-C  bond  with  no  barrier  (Figure  16). 


Minimization  using  previous  MD 
structure  (normal  ReaxFF) 


l=> 

Increase  of  the  Ni- 
C  sigma  bond 
energy  parameter 
by  60% 


e 


C-C  bond  order  decrease 
from  0.91  to  0.71  with  a 
modified  ReaxFF 


Figure  16.  Influence  of  Ni-C  bond  strength  on  C-C  bond  breaking. 


Using  this  modified  ReaxFF  force  field,  we  simulated  new  trajectories  in  which  butane  was 
reacted  with  Ni  at  500K,  600K  and  700K  for  50ps.  In  these  simulations  multiple  CC  cleavage 
events  were  observed,  corroborating  that  the  surface  with  the  strengthened  Ni-C  bonds  would 
decompose  butane  at  room  temperature  (Figure  17). 
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Figure  17.  ReaxFF  MD  simulations  of  butane  conversion  on  the  Ni*(lll)  surface. 


Next,  we  probed  other  metals,  such  as  Pt  and  Ti  which  have  higher  metal-carbon  bond 
strengths.  Using  ReaxFF,  butane  molecules  were  allowed  to  react  with  Pt  and  Ti  surfaces 
(Figures  18  and  19).  As  the  temperature  over  the  Pt  surface  was  increased,  the  first  CC  bond  is 
broken  at  1250K,  similar  to  the  reactions  on  Ni.  For  the  Ti  surface,  CC  cleavage  is  observed  at 
900K.  Therefore,  the  C-C  bond  cleavage  temperature  decreases  by  ~500  K  compared  to  that  for 
the  Ni  surface  and  the  (OH)3Ti-CH3  bond  energy  is  55  kcal/mol  according  to  the  ReaxFF 
potential. 


Figure  18.  Butane  on  the  Pt  (111)  surface,  T=1250-1500  K. 


4.2  Results  &  Discussion 

Despite  some  improvement  in  the  temperature  for  the  catalytic  reactions,  neither  Pt  or  Ti,  are 
is  suitable  for  the  alkane  oxidation  at  low  temperatures. 
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24.4  ps 


24.8  ps 


Figure  19.  Butane  on  the  Ti  surface,  T  =  900  K. 


To  examine  the  oxidation  process  of  the  C-C  cleavage  products  at  low  temperature,  we 
simulated  butane  reactions  on  various  surfaces  (NiO  surface  and  YSZ-Ni  interface)  using  the 
modified  ReaxFF  force  field  (with  increased  Ni-C  bond  strength).  Here  we  used  our  recently 
developed  force-biased  Monte  Carlo  (fbMC)  to  decrease  the  computer  time  to  obtain  reactive 
events.  In  fbMC,  atoms  are  displaced  by  Sr  a  temperature-dependent  component  in  the  force 
direction  for  each  step.  Here  a  controls  the  temperature  dependence  of  the  jump,  kB  is  the 
Boltzmann  constant,  T  is  the  temperature,  Fr  is  the  magnitude  of  the  force,  and  s  is  a  scaling 
factor. 


IkjJ 


a 


leFrl\ 


(1) 


This  gives  a  kick  toward  the  transition  state. 

For  the  NiO(lOO)  surface,  we  observed  C-C  cleavage  and  H20  formation  at  800  K.  We  also 
observed  C-0  bond  formation  (Figure  20). 
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Figure  20.  Butane  on  the  Ni(100)  surface. 

We  increased  the  temperature  of  our  simulations  to  900  K  and  did  ReaxFF  MD  for  ~2  ns  on 
the  NiO  surface.  C-C  cleavage,  H2O  and  C-0  bond  formation  were  observed  at  this  temperature. 
We  still  did  not  observe  formation  of  isolated  CO  or  CO2.  All  the  C  species  remain  bonded  to  the 
Ni  surface.  We  concluded  that  the  Ni-C  bond  is  too  strong  to  fonn  CO  or  CO2.  Thus,  we 
decreased  the  Ni-C  bond  strength  from  65.04  kcal/mol  to  56.8  kcal/mol  and  then  to  48.9 
kcal/mol  and  reran  the  simulations.  The  CO  bond  can  form  on  the  surface  and  desorb  to  the  gas 
phase  when  Ni-C  bond  is  48.9  kcal/mol. 

To  examine  the  CO  oxidation  process,  40  CO  molecules  were  placed  on  the  surface  to 
perfonn  ReaxFF  MD  simulations.  Two  chemical  pathways  to  form  CO2  were  found.  In  the  first 
process,  CO  adsorbs  on  a  surface  Ni  atom,  then  is  oxidized  by  a  neighboring  O  atom  and  desorbs 
into  the  gas  phase.  We  also  observe  CO2  formation  via  the  reaction  of  a  gas  phase  CO  with  an 
under-coordinated  oxygen  bridging  two  Ni  atoms. 

We  also  simulate  reactions  of  the  JP-8  logistic  fuel  on  the  Ni(l  1 1)  surface  (Figure  21).  The 
mixture  representing  JP-8  consists  of  2  isooctane  CgFfig;  4  methyl-cyclohexane  C7H14;  3  m- 
xylene  CsHio;  6  77-dodecane  C12H26;  1  tetralin  C10H12;  4  tetradecane  C14H30  [3].  The  temperature 
was  increased  from  50K  to  2000K  over  68ps.  Similar  to  the  butane  case  C-C  bond  breaking  is 
observed  for  JP-8  molecular  on  Ni  the  (111)  surface.  Both  for  JP-8  and  butane  on  the  Ni(l  1 1) 
surface,  the  main  products  are  H,  FF  from  the  hydrocarbons  reaction  with  Ni  surface. 
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Initial  model  T  =  2000  K,  25  ps 

Figure  21.  ReaxFF  simulation  of  JP-8  on  the  Ni  (111)  surface. 


To  simulate  a  surface  containing  oxygen  that  could  lead  to  C-0  bond  formation  following 
hydrocarbon  decomposition,  the  JP-8  mixture  was  reacted  with  a  surface  containing  both  Ni  and 
yttria-stabilized  zirconia  (YSZ)  (Figure  22).  No  reactions  were  observed  below  2000K  within 
50ps. 

We  simplified  the  system  by  excluding  the  YSZ  electrolyte,  introducing  oxygen  atoms  into 
simulations  of  butane  and  JP-8  decomposition  on  the  Ni  surfaces  and  built  the  models  in  which 
JP-8  constituents  and  Oo  molecules  were  separated  by  a  Ni  slab.  The  parameters  governing  Ni-C 
bonds  were  those  developed  earlier,  in  which  the  Ni-C  bond  energy  was  increased  to  65  kcal/mol 
(the  modified  ReaxFF  with  Ni*  metal)  to  accelerate  hydrocarbon  decomposition.  The  Ni  (411) 
surface  was  used  to  provide  possible  oxygen  diffusion  to  the  surface  after  the  C-C  cleavage 
(Figure  23).  Following  hydrocarbon  decomposition,  oxygen  was  allowed  to  dissociate  on  one 
side  of  the  slab  and  diffuse  to  the  other  side,  where  it  reacted  with  JP-8  decomposition  products. 
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Reflective  waj[ 


Figure  22.  ReaxFF  MD  simulation  of  JP-8 
on  the  YSZ/Ni  interface. 


Reflective  wall 


Figure  23.  JP-8  on  Ni  (411)  surface  +  O2. 


At  600K,  oxygen  atoms  diffused  across  the  Ni  slab  and  reacted  with  surface-bound 
hydrogen  to  produce  water.  Hydrogen  also  diffused  into  the  slab.  The  O  atom  reacts  with  one  H 
to  form  OH  radical  and  then  OH  was  further  oxidized  and  reacted  with  another  absorbed  H  to 
form  H2O.  CH,  CH2,  C2H4  are  formed  via  C-C  cleavage  on  the  Ni  surface.  The  most  abundant 
species  are  CH  and  CH2. 

At  500K,  C-C  cleavage  was  again  observed  on  the  slab.  Oxygen  reduction  to  H20  similar  to 
that  at  600K  was  observed  but  in  smaller  amounts.  Hydrocarbon  species  (CH,  CH2,  CH3)  form 
via  C-C  cleavage  on  the  Ni  surface.  Observation  of  CH3  species  indicates  less  chemistry  than  at 
600  K.  During  the  oxidation  process,  H2O  molecules  and  HO  radicals  formed  on  the  surface  and 
in  the  gas  phase  at  T  =  500  and  600  K.  The  H2O  forms  via  two  oxidization  steps. 

At  400K,  oxygen  did  not  diffuse  through  the  slab,  so  no  hydrocarbon  oxidation  was 
observed. 

4.3  Conclusions 

Summarizing  these  results, 

•  for  low  temperature  activation  of  CC  bonds  we  require  M-C  bond  energies  of  65  kcal/mol 
(similar  to  Ni) 

•  for  low  temperature  activation  of  CH  bonds  we  require  M-H  bond  energies  of  62  kcal/mol 
(similar  to  Ni) 
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•  for  low  temperature  formation  of  CO  and  dislocation  of  C02  bonds  we  require  M-C  bond 

energies  of  49  kcal/mol  (similar  to  Os  and  Ru)  7  and  MO  bond  energies  of  60kcal/mol 

The  conclusion  that  two  different  M-C  bond  emerges  are  required,  one  for  activation  of  the  CC 
bonds  and  the  other  for  fonnation  of  CO  and  C02  products  may  be  achievable.  For  example  in 
the  Ru-Pt  catalysts  for  methanol  fuel  cells  6  the  Pt  atoms  preferentially  break  the  CH  bonds  to 
form  Pt-H  and  Pt-COH  while  the  Ru  atoms  active  H20  to  form  Ru-H  and  RuOH.  The  Ru-OH 
leads  to  formation  of  RuOxHy  clusters  on  the  surface  that  reacts  with  Pt-COH  to  form  C02 
while  the  RuOxHy  conducts  away  the  proton. 

Summarizing,  the  top-down  approach  has  been  used  to  detennine  parameters  needed  to 
permit  low  temperature  alkane  oxidation.  The  next  step  would  be  ReaxFF  simulations  with 
metals  possessing  the  above  bond  energies. 
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5.0  Conference  on  low  temperature  fuel  cells  for  JP8 

We  have  organized  a  conference  on  low  temperature  fuel  cells  with  half  a  day  on  Aug.  9  and  a 
full  day  on  Aug.  12.  Participants  include 

Prof.  Mahdi  Abu-Omar,  Purdue  University,  mabuomar@purdue.edu 
<mailto :  mabuomar@purdue .  edu> 

Prof.  Andy  Herring,  Colorado  School  of  Mines,  aherring@mines.edu 

Prof.  Plamen  Atannasov,  University  of  New  Mexico,  plamen@unm.edu 
<mailto:plamen@unm.edu> 

Prof.  Sanjeev  Mukerjee,  Northeastern  University,  s.mukerjee@neu.edu 

Prof.  Cliff  Kubiak,  University  of  California.  San  Diego,  ckubiak@ucsd.edu 
<mailto:ckubiak@ucsd.edu> 

Prof.  Wenzen  Li,  Michigan  Technological  University  wzli@mtu.edu  <mailto:wzli@mtu.edu> 

Prof.  Andrew  Gewirth,  University  of  Illinois  at  Urbana-Champaign,  agewirth@illinois.edu 
<mailto :  agewirth@illinois .  edu> 

Prof.  William  Tolman,  University  of  Minnesota,  wtolman@umn.edu 

Dr,  James  Gerken  representing  Prof.  Shannon  Stahl,  University  of  Wisconsin, 
stahl@chem. wise . edu  <mailto : stahl@chem. wise .edu> 

Prof.  Ayusman  Sen,  Pennsylvania  State  University 

Manny  Soriaga,  Caltech,  msoriaga@caltech.edu 

Professor  Sri  Nayayan,  Dept,  of  Chemistry,  USC 

Professor  Sossina  Haile,  Dept,  of  Materials  Science,  Caltech 

Dr.  Charles  Hays,  JPL 

Professor  Wenzen  Li.  Michigan  Technical  U,  wzli@mtu.edu 

Professor  Andrew  Herring,  Colorado  School  of  Mines,  aherring@mines.edu 

Professor  Chuan-Jian  Zhong,  SUNY  Binghamton,  cjzhong@binghamton.edu 

Dr.  Hongsen  Wang,  Cornell,  Reprsenting  Prof.  Abruna 

Professor  William  A.  Goddard,  Caltech 

Dr.  Robert  Nielsen,  Caltech 

Dr.  Boris  Merinov,  Caltech 
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